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Abstract 

We manufactured pulsed illuminators emitting in the far infrared 
for the Planck-HFI bolometric instrument ground calibrations. Spe- 
cific measurements have been conducted on these light sources, based 
on Carbon fibers, to understand and predict their properties. We 
present a modelisation of the temperature dependence of the thermal 
conductivity and the calorific capacitance of the fibers. A compari- 
son between simulations and bolometer data is given, that shows the 
coherence of our model. Their small time constants, their stability 
and their emission spectrum pointing in the submm range make these 
illuminators a very usefull tool for calibrating FIR instruments. 

1 Introduction 

We have been studying Carbon fibers as illuminators for FIR instrument 
calibration between 2000 and 2005 with the specific goal of measuring the sum 
of the electrical and optical crosstalk ^ of the Planck-HFI instrument[l][2]. 
Our requirements were threefold: 

• to get a signal ranging up to a few pW detected on the Planck-HFI 
bolometers for the concerned frequencies (ranging from 100 to 857 GHz) 
including transmission efficiency of the cold optics horns and integra- 
tion over the 30% wide frequency bands. 



^The optical crosstalk is caused by potential leakages between the different cryogenic 
stages of the optics facing the bolometers. 
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• to get a relatively small time constant (smaller than 10ms) and a repro- 
ducible signal in a pulsed regime, allowing to stack the measurements 
in order to increase the signal over noise ratio. 

• to be able to "lighten" the fibers without introducing any parasitic 
signal due to EMI-EMC from the electrical pulse used to drive them. 

We did fulfill these requirements. Even more, on the time constant issue, we 
did get better than expected results and the fibers were one of the sources 
used to characterise the time response of the instrument. 

Since the carbon fibers are composite materials, no measurement of the 
particular properties of the ones we used have been made. We therefore focus 
in this paper on their characterisation. In part 2 we give the heat equation 
and the corresponding approximations made in this paper to analyse the 
fiber behaviour. Part 3 describes the experimental setup and in part 4 we 
present the measurements of the relation between resistance and temperature 
and the analysis that leads to the measurement of the thermal conductivity 
and the calorific capacitance. Part 5 describes the fiber modelisation and a 
comparison of simulation results with data. 

2 General framework 

2.1 Equations and assumptions 

The modelisation of the carbon fiber behaviour relies on the one dimensionnal 
heat equation (along the fiber): 



• K is the thermal conductivity, Cp the calorific capacitance, /j, the resis- 
tivity. We assume here that all these parameters do depend on tem- 
perature T. 

• p the mass per unit of volum, and S the surface of a slice [4]. The 
resistance R is obtained through //(T) integration along the length L 
of the fiber. For the fibers we used, L ~ 1mm and the diameter of a 
slice is 6 microns: a photography is shown on figure 1. This picture 




(1) 



where: 
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also shows the way it is thermahsed at each end on a kapton circuit 
through Ag lacquer contacts. 



The fiber is powered by a current (I) generator. 

The radiated power represents less than 1% of thermal losses, even in 
the worst case ^, so it is neglected with respect to conduction. 





Figure 1: Photography of one of the carbon fiber used for the Planck-HFI 
calibration: the diameter is 6fim and the length ^ 1 mm. 



2.2 Basic properties 

The fiber is thermalised at each end at a temperature Tq (cf. figure 1). In the 
steady state regime, and in the simplest case where no parameter exhibits 
any temperature dependence, Eq. (1) gives: 

leading to (the fiber center being the x axis origin): 

and a mean temperature increase: 

^Tmes = J (T{x)~To)dx = ^ . (4) 

L J-L/2 llub 



^Thc fiber radiates about lOfj,W when heated at about 200K by IV pulses, dissipating 
about ImW Joule power. 
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When one stops heating the fiber, the system relaxes to the Tq temper- 
ature, and the transcient regime can be described by the generic Fourier 
expansion solution: 

oo 

T{x, t) = To + Y, AT^e^-^/^"^ cos((2A; + l)7rx/L) , (5) 

fe=0 

where the time constants are given by: 

~ «;(2A; + 1)%2 ' 

the leading term being by far the k — Q component, which corresponds to 
the longest time constant: 

The fiber time constant is then proportional to the square of its length. 
Therefore, smaller time constants may be obtained with shorter fibers, at 
the price of a flux loss. L should stay sizeable with the longest needed 
wavelength, since the flber acts like an antenna. 

For small values of /, AT^es is small, and these formula apply even if k 
varies with T. We use this approximation in the following. 

Resistivity dependence with T can be treated analytically in the simplest 
case where fi;(T) is constant, assuming /^(T) varies linearly as: 

/x(T)=/io[l-a(T-To)] , (8) 

then, Eq. (1) exact solution for the mean temperature increase in steady 
state regime reads [4]: 



AT 



1 

a 



1 



mi \2S \ K 



(9) 



3 Experimental setups 

The measurements described in the following sections have been performed 
within the Saturne cryostat^ providing us with a thermalisation temperature 

^This cryostat at lAS-Orsay has been used for the cahbration of the Planck-HFI in- 
strument. 
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ranging from 300K to 1.7K (we took advantage of the cooling down of the 
setup). Two kind of data are analysed: 

• resistance (or resistance variations) measurements at the edges of the 
fibers done using the dedicated electronics described below for different 
thermalisation temperatures and for small values of the current applied 
on the fibers. 

• bolometers' data translated in terms of incident power for a thermali- 
sation temperature around 2K and a wide range of current values. 

3.1 Dedicated electronics description 




Figure 2: Figure illustrating the electronic circuit used to characterise the 
resistance of the carbon fibers. 

The electronics used for resistance and resistance variation measurements 
is described in figure 2. The signal Vq entering the circuit is a low impedance 
square signal which amplitude, frequency, and width can be changed. The 
output is simply given by: 

V,{t) = {1 - ^^)GVo{t) (10) 

where Rtot corresponds to the resistance of the fiber in series with 15 meters 
long wires. Before any measurement we balance the system in order to get: 
Ri ~ Rtot{t — 0). Taking t — a,t the beginning of the pulse, one gets: 

V,{t) = Vi(0) + ^^(%,er(t) - Rfiberm (11) 
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Figure 3: Example of Vi(t) when a square signal Vq of 0.2 V is applied: on 
the left when r(~ 2.6Q) is switched on and on the right when it is off. The 
device reaches a precision of the order of .05f2 on the resistance variation 
measurements. 



The gain factor is measured by switching on and off the small additionnal 
resistor r. 

An example of the signal measured at the edges of the fibers is illustrated 
on figure 3 (cf. section 4.2.1 for more details). 



3.2 Saturne setup 

The setup of the fibers within the Saturne cryostat did take into account 
the constraints of all the measurements which were needed to be done to 
calibrate the Planck-HFI instrument (and which are not discussed here [3]). 
The fibers were therefore installed on a side of a 3 positions wheel: for one 
of them, the fibers were facing the entrance of part of the cold optics of 
the instrument as shown on figure 4. Two additionnal fibers were installed 
behind a hole in a mirror facing the focal plane and illuminating all the 
bolometers synchroneously, for time constant measurements. Being further 
from the instrument, the mirror fibers illuminate HFI horns about 300 times 
less than the wheel fibers. 



6 



Figure 4: Schematics and picture of the setup used for the characterisation of 
the Planck-HFI instrument (which cold optics is represented in yellow on the 
left and which is in gold on the picture) with the carbon fibers illuminators 
(resp. in blue and in grey). 

4 Characterisation 

4.1 Resistance and temperature 

At room temperature the measured resistivity of the fibers is 1.8 10~^Q.m 
which gives for 1mm fibers a resistance of 640 Ohms at 300i^. This resistance 
does depend on the thermalisation temperature as shown on figure 5 for seven 
fibers (on the left): the dispersion of these measurements includes the spread 
of the resistance of the Ag-lacquer contacts and the fact that the length of 
the fibers lies between 0.94 and 1.06 mm. The figure on the right shows, for 
the same seven fibers, their resistance corrected for their value measured at 
80K and divided by the slope ^ at T = 80K to correct for length effects. It 
illustrates the homogeneity of the results. The contact resistance due to the 
Carbon- Ag lacquer interfaces is of the order of lOOQ. 

4.2 Thermal conductivity 

We measure and check the thermal conductivity k, using the two experimental 
setups described in section 3. 
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Figure 5: On the left: resistance of the fibers as a function of the ther- 
mahsation temperature (a polynomial fit is superimposed). On the right: 
comparison of the same data when properly normalized (cf. text). 



4.2.1 With the dedicated electronics 

The first step of the analysis using the dedicated electronics is to estimate the 
resistance variation of the fibers as a function of the thermalisation tempera- 
ture for very small voltages applied on their edges (typically under lOOmV), 
heating the fiber about lOK above Tq. For such a voltage V^, we measure 
the amplitude Vf of the signal between t = and its asymptotic value, and 
the difference AVr of the voltage signals measured at the edges of the same 
fiber when pulsed in series with the resistance r = 2.65Q and without r (cf. 
Figure 3). The variation of the resistance induced by a tension Va applied 
on the fibers is given by: 

AK = ^ r (12) 

We deduce the value of the mean temperature to which the fiber is heated 
inverting the polynomial function used to fit the R(T) behaviour of the fibers 
presented in section 4.1. For instance, above 40K, the resistance is propor- 
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tional to the temperature (cf. section 4.1) and one gets from (4) and (8): 

Va^La 



USARa 



(13) 



The results are shown on figure 6 for seven superimposed fibers, k is well 
described in this temperature range by the parametrisation: 



K = Ko + KaT^ + KgT^ (14) 

the full line curve corresponds to a fit of data with: 
kq ~ 0.3 ± O.lW/Km, K2 ~ (2.5 ± 0.8)10-^W/kV, and /tg ~ (-5. ± 
3.) 10"^ W/kV. 

These results are in good agreement with other measurements on carbon 
fibers[5]. 

The small slope of k(T) at low temperature makes the fiber signal very 
stable even if Tq is not stabilized. For instance, with Tq around AK with 
variation of the order of lOOmi^', the fiber signal varies by less than 10~^. 




Figure 6: Thermal conductivity function of the temperature for seven 

Carbon fibers. The full line corresponds to a fit according to eq. (14) 
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Figure 7: Incident power measured by a lOOGHz Planck-HFI bolometer (in 
pW) - without corrections for optical efficiency of the cold optics of the 
instrument - as a function of the voltage applied at the edges of the fibers for 
2 different fibers: a 1mm length fiber (in black) and a 0.6mm one (in red). 
The superimposed fitted function corresponds to the paramaterization given 
in equation (17). 



4.2.2 Bolometers data 

In order to cross-check this result, one makes use of the Planck HFI instru- 
ment. We pulse the fibers with square signals at approximately 1 Hz, and 
compute the amplitude Sboio of the resulting synchroneous signal measured 
on the bolometers (translated in terms of incident power on the detectors). 
Shoio is proportional to the integral of the temperature increase over the 
length of the fiber ^. Table 1 gives the measured values of the flux (in pW) 
with the bolometers for the fibers installed on the mirror and pulsed with IV 
amplitude signals^. 

Figure 7 shows the power measured by one of the bolometers (at lOOGHz) 
for two "mirror" fibers (a small 0.6mm fiber and a longer 1mm one) as a 
function of the voltage applied on the fiber. 

^One assumes that, being in the Rayleigh Jeans domain, the emission spectrum of the 
fiber is proportional to T. 

^These values have to be multiplied by a factor 350 if one wants to estimate the flux 
when the fiber is directly facing the entrance of the horn of the cold optics in front of the 
bolometers. 
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lOOGHz 


143 GHz 


217 GHz 


353GHz 


545 GHz 


Flux(10-^ pW) 


11. 


10. 


21. 


117. 


27. 



Table 1: Measured flux for each Planck-HFI frequency band for the flbers 
installed on the mirror and pulsed with IV amplitude signals. 



These results can be understood in the hght of the framework given in 
section 2. In the permanent regime, using (14) to parametrize k{T) and 
neglecting R{T) dependence, one gets from equation (1): 

(^T+^^^^:^] ^ (15) 



dx^ \" 3 4 y LSR 
Since T = Tq for x = ±L/2, one obtains after integration: 



KoT{x) + ^T%x) + ^T\x) = A{x)+j (16) 



where: 



A{x) = ^(^V4 - x') , and 7 = ^^o^o + |To=^ + ^To^ . 

Since here Tq ^ 2K, it can be neglected with respect to T. Hence, 7 = 0. 

A{x) simplifies to A{x) = V^/{8n){l - {2x/Ly). Hence, the solution 
T[x) is only a function of V and the reduced variable ^ = 2x/L. 

For smal V values, the first term dominates the left-hand side of (16), 
that reduces to (3), and T{x) has a parabolic shape. At medium V values, 
T{x) solution on most part of the fiber imphes mainly the term, and T{x) 
has a fiatter profile with a -y^l — ^ shape. 

At high V value, the negative Ks term plays a higher role and is responsible 
for the increase in the T{x) slope. 

The phenomenological expression: 

S,oio oc GoV' + ^V'/' + ^V'/' (17) 

is found to fit accurately both experimental (see figure 7) and simulated 
data (see section 5). 

Here, the applied voltages are higher that the ones used in the previous 
section, still the model fits accurately the data and both results are in very 
good agreement. 
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4.3 Time constants and Calorific capacitance 

With the dedicated electronics of section (3.1), and applying only small volt- 
ages to the fibers, one can extract the fiber time constant by fitting the 
exponential behaviour of Vi{t) in a wide range of thermalisation tempera- 
ture Tq. The results on r are illustrated on figure 8: on the left hand side the 
data are in ms while on the right there are normalized to the mean value ob- 
tained for each fiber with T > 250K, showing that this behaviour is coherent 
from one fiber to the others. 
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Figure 8: On the left: Time constant (in ms) as a function of the thermal- 
isation temperature. On the right: same data but normalized to the mean 
value obtained for each fiber with T > 2507^. 

In the temperature range we explore, the fiber time constants present 
a minimum around lOOii', and rise slowly up to 300K. The higher time 
constants measured for a temperature smaller than 50K are responsible for 
long decay times at the end of the pulse. During the Planck HFI calibration 
[3], this drawback has been corrected for using a small permanent current on 
the fiber in addition to the pulse. This permanent current allows to maintain 
most of the fiber at temperatures where the time constant remains small, and 
we recover a fiber decay time in the same range as the rise time. 

Cp can be deduced from the measurement of the time constant and from 
the K modelisation through relation (7). Figure 9 shows individual results 
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Figure 9: Calorific capacitance Cp as a function of the temperature for several 
Carbon fibers. 
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for seven fiber calorific capacitance. Parametric adjustments of the standard 
type Cp/T = Cpo + T^Cp2 are superimposed. A good qualitative agreement 
is met with expected values for Carbon fiber at low temperature [6]. 

5 Simulation 

We developped a simple simulation code solving numerically equation 1. Us- 
ing the Crank Nicholson scheme [7] and introducing the fitted dependence 
on the temperature of R, k and Cp, we can compute the behaviour of the 
mean temperature of the fiber as a function of the input voltage: it is shown 
on figure 10 which reproduces the results measured on the bolometers (cf. 
Figure 7). The parameters extracted from data with some approximations 
are therefore coherent with the simulations. 
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Figure 10: Simulation of the mean temperature of the fiber as a function of 
the amplitude of the applied voltage for several fibers lengths. 
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6 Conclusion 



This article describes the physical parameters of the carbon fiber illuminators 
used in the calibration of the Planck-HFI instrument. We have shown that, 
thanks to their low time constant ( < 10ms for 1mm fibers at 2K), they can 
be used in a pulsed regime without introducing any electrical parasitic signal 
on HFl's bolometers, and that their emission spectrum gives a significant 
amount of signal in the submm domain. 

We have detailed an analysis of the temperature dependence of the resis- 
tance, the calorific capacitance Cp, and the thermal conductivity «, of theses 
fibers based on measurements from 300 to 1.7K. Wc end up showing the 
good agreement between simulations of the ID heat equations making use of 
the extracted parameters (-R(T), Cp{T) and n{T)) and data measured of the 
lOOmK HFI's bolometers. 

This well understood picture of the fibers behaviour makes them a very 
useful! tool for FIR instrument calibration. 
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